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ABSTRACT;  A  method  of  estimating  source  energy  levels  for  shallow  underwater 
explosions  used  as  acoustic  source;  la  developed.  The  method  uoes  Idealized 
quail-theoretical  preeaure-tlce  recorde  and  digital  process!.-.,;  techniques.  The 
o«e<l  for  such  a  aethod  of  estimating  source  levels  Is  discussed  in  light  of  the 
fact,  that  when  the  direct  arrival  pressure  pulse  is  not  tlno-separable  from  the 
aurface-reflaoted  arrival,  no  direct  measurement  of  source  level  is  pocaible. 

Source  levels  are  calculated  for  )  lb  of  TNT  at  charge  depths  of  60  and  5C0  feat 
and  1.8  lb  of  TNT  at  depths  of  60,  ;00,  and  £C0  feet.  The  1.8-lb  results  epcly 
to  the  widely-used  underwater  sound  signals  (uNS)  K<  57,  Kk  61,  and  hit  32  with 
appropriate  nominal  depth  settings;  the  3-lb  charge  results  correspond  to 
configuration*  used  in  an  extensive  research  study. 

The  Influence  of  depth  variation,  surface  reflections,  and  resolution  bandwidth 
has  beos  examined.  It  Is  shown  that  seemingly  small  variations  in  burst  deptn 
can  cause  variations  up  to  6.3  do  for  certain  narrow-band  3ourco  levels.  The 
effect  of  the  surface  reflection,  sometimes  erroneously  included  as  part  cf  the 
aource  level,  is  demonstrated. 

Comparisons  of  the  Idealised  model  ef  this  report  with  Weston's  analytical  model, 
the  experimental  results  of  Christian,  ar.d  the  experimental  results  of  Turner  and 
■SoriEger  show  maximum  differences  of  8,  3,  and  5  db,  respectively,  in  6o=e 
instances. 

All  discrepancies  have  not  end  cannot  be  resolved  vlthout  further  conparl.-cr.s  with 
high-quality  broadband  ntaeurecents  cf  shallow  underwater  explosions.  It  Is 
tentatively  recommended  that  shallow  source  levels,  where  the  surface  reflection 
cannot  bo  separated  out,  be  determined  by  the  method  of  the  present  study. 
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Prepared  byi 

Joel  8.  Gaapin  and  Verna  K.  Shuler 


ABSTRACTS  A  method  of  estimating  eource  energy  levels  for  shallow  underwater 
explosions  used  as  acoustic  sources  la  developed.  The  aathod  uses  idealised 
quasi-theoretical  presaure-tiae  records  and  digital  processing  techniques*  The 
need  for  aueh  a  method  of  estimating  source  levels  is  discussed  in  light  of  the 
fact  that  when  the  direct  arrival  pressure  pulse  is  not  time-separable  from  the 
surface-reflaoted  arrival,  no  direct  measurement  of  source  level  is  possible. 

Source  levels  are  calculated  for  3  lb  of  TNT  at  charge  depths  of  60  and  50 0  feet 
and  1,8  lb  of  TNT  at  depths  of  60,  300,  and  800  feet.  The  1.8-lb  results  apply 
to  the  widely-used  underwater  sound  signals  (SUS)  Mk  57,  Mk  61,  and  Mk  82  with 
appropriate  nominal  depth  settings;  ths  3- lb  charge  results  correspond  to 
configurations  used  in  en  extensive  reeearoh  study. 

The  influence  of  depth  variation,  surface  reflection!,  and  raaolutlon  bandwidth 
has  baen  axamlned.  It  is  shown  that  ssamlagly  maall  variations  in  burst  dspth 
can  cause  variations  up  to  6.3  db  for  oartain  narrow-band  source  levels.  The 
affect  of  tha  surface  reflection,  sometimes  erroneously  included  aa  part  of  the 
source  level,  la  demonstrated. 

Comparisons  of  the  idealised  model  of  this  report  with  Weston's  analytical  model, 
the  experimental  results  of  Christian,  and  the  experimantal  results  of  Turner  and 
Scrimger  show  maxlmm  dlffarensts  of  8,  3,  and  5  db,  respeetlvsly,  in  some 
instances. 

All  discrepancies  have  not  and  oennot  be  resolved  without  further  comparisons  with 
high-quality  broadband  measurements  of  shallow  underwater  explosions.  It  is 
tentativsly  reooamended  that  shallow  souroa  levels,  where  the  surface  reflection 
cannot  be  separated  out,  be  determined  by  the  method  of  the  present  study. 
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This  report  is  part  of  a  continuing  study  of  the  near-field  acoustic  characteristics 
of  underwater  explosions.  Snail  explosive  charges  are  frequently  the  sources  used 
in  underwater  acoustic  experiments,  particularly  those  deeding  with  long-range 
propagation  of  low-frequency  sound.  Such  acoustic  experiments  rarely  Include 
direct  near-field  measurement  of  the  charge  output  for  each  shot;  source  level 
values  are  typically  estimated  quantities.  Since  there  is  no  single,  generally 
accepted  method  of  estimating  source  levels  at  present,  different  users  of 
explosion  sources  sometimes  arrive  at  significantly  different  values  of  thlB 
important  parameter. 

This  report  presents  idealized  acoustic  spectra  for  several  explosive  charge 
weight /depth  configurations  frequently  used  in  underwater  acoustics  research 
studies.  These  spectra,  ^ftilch  were  derived  from  quasl-theoretical  pressure-time 
signatures,  have  not  yet  been  adequately  validated.  Additional  work  required  to 
confirm  or  modify  the  Idealized  source  level  functions  is  outlined. 

This  work  was  performed  on  ONR  Research  Project  R  2406  (Project  Order  1-0062). 

The  authors  are  Indebted  to  Miss  Ermine  A.  Christian  for  her  technical  guidance  and 
discussions  of  this  work. 


ROBERT  WILLIAMSON.  II 
Captain,  U5N 
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SOLACE  LEV'S  0?  3HALI/G*  UNUKRWATER  EXPTiXIONS 


1.  INTRODUCTION 


1.1  Because  of  the  widespread  use  of  underwater  explosions  as  acoustic  sources 
for  long-range  propagation  studies,  tbe  problem  of  the  determination  of  explosion 
source  energy  levels  for  underwater  explosions  has  became  increasingly  important  in 
recent  years.  Transmission  loss  Is  perhaps  the  foremost  problem  In  underwater 
acoustics  research,  and  reliable  source  levels  are  needed  if  transmission  loss  is  to 
be  found  from  measurements  of  received  signal  level.  For  many  explosion  geometries, 
a  more  or  less  direct  measurement  of  the  source  level  Is  possible.  When  the 
explosive  source  Is  so  shallow  that  the  perturbations  due  to  the  water  surface  sue 
not  time- separable  from  the  direct  effects,  no  direct  measurement  of  the  source 
level  can  be  made.  An  analytical  approach,  using  quasi -theoretical,  Idealized 
oressure-tlme  curves  to  determine  the  energy  spectrum,  can  be  useful  for  calculating 
the  source  levels  of  shallow  explosions.  One  such  calculation  1,  the  standard 

tr  ecu'v.t  Weston, 

l.?  The  work  documented  in  this  report  uses  idealized  pressure-til^  curves 
Much  comprise  a  more  detailed  representation  of  the  free  water  curves  than  those 
used  by  Weston.  Using  a  digital  computer,  and  fast  Fourier  transform  techniques, 
the  energy  spectra  are  determined,  and  may  be  Integrated  in  any  desired  analyslr 
bands.  This  technique  has  been  used  to  determine  the  source  levels  at  a  100-yard 
reforence  range**  for  several  conditions  of  interest  in  order  to  demonstrate  the 
method  as  well  as  to  provide  source  levels  for  these  conditions.  These  source 
levels  are  displayed  graphically  for  frequencies  up  to  500  Hz  and  tabulated  In 
octave  and  l/3-octave  bands  for  frequencies  of  11  Hz  to  11  kHz, 

1.3  The  discusalon  below  Includes  a  brief  treatment  of  underwater  explosion 
phenomenology  related  to  the  vise  of  explosions  as  acoustic  sources,  as  well  as  a 
discussion  of  the  problems  encountered  in  the  use  of  explosive  sources.  Pertinent 
previous  work  on  explosive  source  levels  is  reviewed  and  related  to  the  present 
calculations.  The  details  of  the  construction  of  the  idealised  pressure-time 
curves  and  the  results  of  the  analysis  are  given.  Some  preliminary  comparisons  of 
the  present  results  with  those  of  previously  used  techniques  have  been  made  and 
evaluated 

2.  UNDERWATER  EXPLOSION  PHENOMENOLOGY 

2.1  M&ny  of  the  problems  encountered  In  dealing  with  underwater  explosions  as 
acoustic  sources  are  due  to  the  complex  nature  of  the  explosion  phenomena,  and 
consequently,  the  pressure-time  signal.  The  processes  Involved  are  fairly  well 
understood.  The  following  brief  discussion  vlll  touch  only  on  the  points  relevant 
to  the  pressure  signal  produced.  The  reader  is  referred  to  the  book  by  Cole^  and 
co  Weston’ j  pa perl  a9  good  general  references. 

Refers  to  references  on  page  2k  and  25. 

**  The  100-yard  reference  range  Is  chosen  following  Weston  (Ref.  1),  Christian 

(Ref.  2),  end  others.  The  scaling  of  energy  spectra  from  one  range  to  another  is 
discussed  by  Christian  in  Ref.  3. 
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2.2  At  the  time  of  detonation  of  an  underwater  exploelve  charge,  the  explosive 

la  transformed  tor  the  detonation  mere  into  a  nail  sphere  of  gaseous  explosion 
products  at  high  pressure,  and  a  shock  wv«  is  propagated  out  sard.  The  gas  globe, 
being  at  a  vtimrt  nuk  higher  ♦-*»*»  the  lurraundlne  Milte.  heelne  in 

site.  As  the  sphere  emends,  the  pressure  in  the  tail  at  the  shook  wave  Is  reduced, 
corresponding  to  ths  decrease  in  preseure  in  the  bubble.  Tbo  bubble  oontinues  to 
expand  until  it  reaches  its  maximue  extent — »t  >hld>  time  the  traasnitted  pressure, 
herring  become  negative  with  respect  to  shbient  hydrostatic  pressure,  is  at  a 
minimum.  The  hydrostatic  pressure  Is  now  greater  than  that  In  the  bubble; 
consequently,  the  bubble  begins  to  coot  rant  and  Is  carried  by  momsntun  through  the 
equilibrium  vo liras  to  a  very  snail  minima  A  positive  pressure  pulse,  the  first 
bubble  pulse,  le  emitted  end  the  bubble  expansion  and  contraction  cycle  begins 
again  and  nay  continue  through  as  many  as  ten  cycles.  Ieoh  succeeding  cycle  takes 

a  shorter  time,  and  each  succeeding  bubble  pulse  has  a  loner  peak  pressure  and  less 
energy.  Figure  1  shone  a  typical  explosion  pressure  signal  along  with  a  depiction 
of  the  bubble  site  at  various  stages  at  ths  cycle. 

2.3  The  tins  interval  between  the  shock  nave  end  first  bubble  pulse,  termed 
the  first  bubble  period,  Is  a  sensitive  function  of  charge  weight  and  depth.  If 
the  effect  of  the  boundaries  of  ths  water  Is  neglected,  the  first  bubble  period, 

T^,  is  given  by  the  expression 

Tt  -  KW1//3/Z5//6|  (1) 

idiere  W  is  the  charge  weight,  Z  is  the  total  hydrostatic  bead,  and  K  is  a 
coefficient  depending  on  the  explosive  composition.*  Some  representative  values 
of  for  given  conditions  are  ehown  In  Table  1  (see  also  Figure  2). 

Table  1 


First  Bubble  Period,  T^,  for  TNT  (seconds) 


Charge  Weight 

anrir.  ini 

_ life! _ 

22 

222 

2220 

1 

.11 

.024 

.0036 

10 

.24 

.050 

.0078 

100 

.51 

.11 

.017 

The  functional  dependence 

of  the  various 

explosion  pressure  history 

parameters  is  discussed  in  section  5,  below. 

2.5  While  the  explosion  bubble  is  pulsating,  It  aay  experience  significant 
upward  migration  as  a  result  of  buoyancy.  Thu,  ths  explosion  heoomss  a  moving 
source,  with  succeeding  bubble  pulses  originating  at  shallowtr  depths.  Bubble 
migration  has  been  treated  by  8nay3 .  Calculations  show  that  for  the  charge  sizes 
and  geometries  considered  in  this  report,  migration  is  not  a  significant  effect. 

For  the  3-lb,  60-ft  condition,  the  strongest  migrating  condition  considered  in  this 
report,  the  bubble  migrates  about  2  feet  to  the  2nd  bubble  maximum.  This  is  not 
considered  s  mijar  influence  on  the  source  levels.  Far  larger  charges,  at  modarmts 
depths,  the  effect  is  more  important. 


may 


exhibit  a  slight  dependence  on  charge  depth  tar  some  explosives. 
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2.6  The  explosion  pressure  signal  is  in  general  modified  by  the  effect  of 
the  surface  of  the  water.  The  reflected  signal  generally  exhibits  \  phase 
change  of  180°  vlth  respect  to  the  direct  arrival.  For  shallow  explosions, 
this  reflected  signal  is  not  time  separable  frot  the  direct  arrival-  in?*  relation- 
Cwip  0^  toe  surface  reflected  pulse  to  the  study  of  explosive  source  levels  is 
discussed  in  the  next  section. 

3.  UNDERWATER  EXPL06JQIC  AS  ACOUSTIC  SOURCES 


3>1  An  underwater  explosion  is  e  broadband  energy  source,  with  considerable 
acoustic  energy  radiated  over  a  broad  range  of  frequencies.  Explosives  are 
generally  inexpensive  and  produce  a  great  smount  of  energy  per  dollar  in  any 
frequency  band  of  Interest.  Underwater  explosions  may  be  easily  detonated 
at  any  desired  depth,  without  the  need  far  connecting  vires  or  cables.  Explosive 
sources  are  non  directional;  they  radiate  a  large  amount  of  energy  in  all 
directions.  This  is  both  an  advantage  and  a  drawback  of  explosive  sources. 

The  explosion  may  radiate  adequate  energy  in  the  desired  direction  for  a  particular 
application,  but  the  fact  that  it  radiates  the  seas  energy  in  all  directions 
leads  to  hlgft  reverberation  noise  levels.  Another  difficulty  Is  that  the  explosive 
source  is  not  well  suited  to  get  repeated  "looks"  at  the  target  in  active  sonar 
applications.  In  spite  of  these  drawbacks,  the  explosive  source  has  found  wide 
usage  in  the  acoustics  community.  The  general  subject  of  underwater  explosions 
as  acoustic  sources  is  discussed  by  Vsstorr-  and  Uriek. 

3>2  The  source  level  of  an  underwater  explosive  charge  gives  a  measure  of 
the  acoustic  energy  in  a  specified  frequency  band  that  is  radiated  by  the  charge 
to  a  particular  point  in  spaces  it  is  a  frequency-domain  representation  of  the 
pres sure- time  function  shown  in  figure  1.  In  many  cases,  experimental  conditions 
are  such  that  the  boundaries  or  the  medium  itself  Introduce  perturbations  in 
the  pressure  signal  before  the  bubble  pulses  have  bubslded.  Here  we  are  concerned 
with  such  perturbations  as  those  caused  by  the  phase-shifted  surface  reflection. 
When  the  surface  reflection  is  superimposed  on  the  direct  arrival,  neither  direct 
calculations  of  the  source  level  from  the  measured  pressure  signal  nor  analogue 
measurement  of  the  source  level  Is  possible.  Because  the  surface  reflection  Is 
not  time-separable  from  the  direct  effects  of  the  explosion,  same  manipulation, 
such  as  the  frequency  demodulation  of  Parkins, '  is  necessary  bafore  a  true  source 
level  can  be  determined  from  recordings  in  iAicn  the  rarefaction  wave  interferes 
vlth  the  direct  wave.  The  time-domain  treatment  discussed  here  is  another  approach 
to  the  problem  of  accurately  determining  source  levels  of  "shallow"  explosions. 

3.3  In  the  context  of  this  report,  "deep"  explosions  are  those  for  >&lch 
the  surface  reflection  arrives  after  the  bubble  pulees  have  subsided.  For  these 
"deep"  explosive  sources,  there  is  no  confusion  as  to  that  is  meant  by  source 
level.  It  involves  the  total  acoustic  energy  radiated  by  the  charge  to  some 
fixed  range  (euch  as  IOC  yards).  For  "shallow"  sources,  this  concept  remains 
valid;  however,  any  measurements  of  the  source  levels  of  such  "shallow"  explosions 
unavoidably  contain  energy  contributions  due  to  the  presence  of  the  surface. 

Some  workers  have  erroneously  modified  the  above  concept  of  source  level  to 
Include  surface  reflection  effects  since  no  direct  determination  of  "shallow” 
source  levels  can  otherwise  be  made. 
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I.1 2*  The  time  interval  between  arrival ■  of  the  direct  end  surface  reflected 
wave  Is  &  geometric  function  independent  of  the  charge  weight .  The  bubble 
periods,  and  thus  the  duration  of  the  direct  arrival,  depend  upou  the  weight 
as  well  as  the  depth  of  the  charge.  It  le  the  relation  akin  «?  the  ti ss  of 
arrival  of  the  surface  reflection  to  the  duration  of  the  direct  pulse  that 
determines  whether  source  level  measurements  can  be  made:  when  the  surface 
reflection  occurs  after  the  last  detectable  bubble  period,  direct  source  level 
measurements  are  possible. 

3.5  Both  60- foot  conditions  of  this  report  are  those  for  which  no  direct 
source  level  measurements  can  be  made,  for  the  other,  deeper  depths  (300,  500 
and  3 00  feet;,  the  observation  point  dictates  the  possibllty  of  soiree  level 
measurements.  Under  many  experimented  conditions,  ons  nay  not  have  any  latitude 
in  selecting  the  observation  point  (such  as  *»en  a  fixed  hydrophone  system  is 
being  utilized). 

3.6  An  analytical  method  to  determine  "shallow"  explosion  source  levels 
has  been  evolved  and  is  the  subject  of  sections  5*7  of  this  report.  The  method 
is  also  applied  to  several  cases  of  "deep"  explosions  to  provide  source  levels 
for  cases  In  which  no  appropriate  near-field  measurements  have  been  made. 

4.  PREVIOUS  WORK  ON  SOURCE  LEVELS 

li.l  The  previous  work  on  the  close-ln  energy  spectra  of  underwater  explosions 
may  be  loosely  classified  into  the  two  following  categories: 

1.  "nalytlcal  treatments  based  on  the  Pourler  analysis  of  simple 
idealized  pressure  waveforms  (usually  combinations  of  exponentials)  and 

2.  Analysis  of  Experimental  Data. 

These  tvo  types  of  treatment  will  be  discussed  briefly  below  end  examples  provided. 

4.2  Analysis  of  Simple  Waveforms. 

4.2.1  These  treatments  take  advantage  of  the  feet  that  the  shock  wave 
from  an  underwater  explosion  has  an  Initially  exponential  decay.  The  shock  wave 
pressure-time  history  is  treated  as  an  exponential  pules, 

F(t)  =■  Poe  (3) 

where  P  ,  the  peak  shock  wave  pressure,  and  ft,  the  decay  constant,  are  usually 
evaluated  from  published  empirical  relationships.  Treatments  using  only  the 
shock  wave  in  this  manner  are  highly  unrealistic  as  the  energy  contribution  of  the 
bubble  pulses  Is  very  important.  Therefore,  bubble  pulses  are  usually  taken 
into  account.  Typically,  the  bubble  jwlse  is  idealized  as  a  double  exponential. 

The  time  of  the  peak  bubble  pressure  is  taken  from  experiment.  Such  simple 
idealized  waveforms  have  the  advantage  of  being  amenable  to  solution  for  the 
energy  spectrum  In  closed  form.  This  allows  insight  into  the  variation  of 
different  spectral  parameters  with  changes  in  the  pressure-time  parameters. 
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k.2.2  The  major  diaad* antage  of  these  treatments  is  that  exponentials 
do  not  accurately  model  the  explosion  pressure  pulse.  While  an  actual 
explosion  shock  wive  has  an  initially  exponential  decay,  after  a  time  nt 
roughly  2m,  the  decay  rate  lessens  markedly.  The  shock  ware  pressure  becomes 
negative  (relative  to  hydrostatic  pressure)  at  a  time  strongly  dependent  on  the 
depth  of  burst  and  contains  a  significant  amount  of  energy  in  the  negative 
phase.  The  bubble  pulses,  having  negative  pressure  phases  both  before  and 
after  the  bubble  maxima,  are  also  not  veil  modeled  by  exponentials.  Nonetheless, 
such  treatments  have  proven  both  practical  and  convenient  for  determination  of 
source  levels. 

4.2.3  The. prime  example  of  this  type  of  work  is  the  vell-knovn 
treatment  by  Weston.*  It  should  be  noted  that  in  his  paper  Weston  presents 
tvo  methods  of  determining  source  levels,  an  analytical  method  and  one  based  on 
experimental  data.  He  recommends  the  use  of  his  experimental  results  as  the 
more  accurate  method  of  obtaining  source  levels.  It  is  apparent,  however, 
from  the  nominal  explosive  conditions  shown,  that  surface  reflections  were 
Included  in  some  of  his  data.  Por  this  reason,  the  following  discussion  will 
consider  only  the  analytical  formulation.*  Weston  considers  a  pressure 
history  consisting  of  a  shock  wave  and  tvo  bubble  pulses.  He  chose  to 
separate  his  analysis  into  tvo  parts,  the  low  and  the  high  frequencies.  At 
low  frequencies,  the  shock  wave  and  bubble  pulses  are  replaced  by  their 
respective  impulses.  In  order  to  correct  the  residual  impulse  of  his  model, 
Weston  introduces  a  negative  Impulse,  corresponding  to  a  steady  negative 
pressure,  which  causes  the  total  residual  impulse  to  be  zero.  At  higher 
frequencies,  Weston  shifts  to  an  exponential  shock  wave  followed  by  tvo 
double  exponential  bubble  pulses.  He  uses  a  formula  for  summing  the  spectral 
contributions  of  the  three  pulses  which  neglects  the  interference  patterns 
between  the  pulses.  Weston's  analysis  Includes  wave  form  variations  only 
through  the  bubble  periods,  and  takes  no  account  of  the  variation  in  energy 
and  impulse  in  the  first  positive  phase  (shock  wave)  with  source  depth. 

Weston's  analysis  is  discussed  further  in  section  7. 

8 

4.2.4  Other  examples  of  this  type  of  analysis  are  given  by  BrurJbach 
and  Sllfko?.  This  latter  work  treats  only  the  shock  wave  portion  of  the 
pressure  time  history. 

4.3  Analysis  of  Experimental  Data 

4.3.1  The  determination  of  source  levels  by  the  direct  analysis 
of  experimental  pressure-time  curves  is  a  widely  used  procedure.  The 
measured  pressure  histories  are  analyzed,  either  with  analog  or  digital 
processing,  to  yield,  more  or  less  directly,  the  energy  spectrum  at  a 
reference  range.  This  method,  while  the  most  direct,  has  several  drawbacks. 

The  measurements  are  rarely  made  at  the  chosen  reference  range.  Instead, 
various  propagation  models  are  used  to  extrapolate  back  to  this  range. 

Models  varying  from  spherical  spreading  through  refractive  calculations  and 
sophist! catec'  finite  amplitude  analyses  have  been  used.  In  addition,  the 
desired  source  level  is  not  generally  identical  with  the  energy  spectrum  of 
an  experimental  pressure  history,  even  if  recorded  at  the  appropriate  range. 


*  When  other  authors  use  Weston's  source  levels,  it  i e  often  not  made  clear 
which  of  his  formulations  has  been  used. 
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The  pressure  history  is  subject  to  distortions,  such  as  reflections  fro®  the 
bat  too  and  surface  of  the  eater  Wiich  are  not  pert  of  the  source  level. 

In  addition,  localised  effects  such  as  bulk  cavitation  nay  have  a  profound 
effect  on  a  eeaeured  spent rue,  eh lie  being  relatively  unimportant  at  other 
r ‘uruMTsort,  toe  accurate  racoramg  of  the  pressure  history 
ceui.es  serious  experimental  difficulties.  High-quality  broadband  (~  0  -  20KJTC) 
recordings  are  difficult  to  achieve.  Many  experimenters  record,  therefore, 
only  In  limited  frequency  bands.  This  makes  it  very  difficult  to  Isolate 
the  spurious  effects  mentioned  above,  with  the  consequence  that  soss* 
experimenters  seem  to  hare  included  such  effects  in  their  source  Levels. 

4.3.2  A  valuable  effort  in  this  direction  le  the  paper  by  Christian2, 
nils  paper  documents  the  measured  spectra  of  a  large  number  of  charges.  The 
recordings  vue  high  quality,  broadband;  the  analysis  mas  digitally  performed. 
The  results  were  filtered  Into  octave  bands  and  presented  in  a  scaled 

form  as  a  family  of  curves  representing  all  charge  weights  at  source  depths 
from  300*  to  ?2, 000  feet.  By  reducing  the  energy  by  the  4/3  power  of  the 
charge  weight,  Christian  collapsed  the  data  from  various  charge  sites  at 
the  same  depth  into  a  single  curve.  The  recordings  were  carefully  analysed 
to  eliminate  spurious  effects.  Since  the  recordings  of  the  data  was  done 
directly  above  the  charge,  refractive  effects  were  eliminated.  Theoretical 
considerations  were  used  to  separate  the  effects  of  charge  depth  and  range 
to  the  measuring  hydrophone.  Christian's  curves  have  been  widely  used  to 
provide  source  Levels  in  situations  where  no  close-in  measurements  have 
been  made,  with  several  workers  tending  to  substantiate  her  results,  notably 
fClbblethlte  and  Denham  ,  On  the  other  hand.  Turner  and  Scrlmgerl2  report 
close-in  measurements  that  do  not  agree  with  Christian's  curves  for  same 
frequencies.  This  point  is  discussed  further  in  section  7.4. 

4.3.3  Other  work  along  these  lines  has  often  concerned  itself  with 
the  determination  of  the  source  spectrum  for  a  condition  of  particular  ^ 
Interest  to  the  marker,  such  as  Maples  and  Ihcrpl3.  Christian  and  Blaik 
give  a  good  general  treatment  based  on  experimental  measurmneBts. 

4.3.4  While  this  type  of  work  has  proven  generally  useful,  the 
direct  determination  of  source  levels  for  very  shallow  explosions  is  not 
possible, as  pointed  out  in  section  3-  Efforts  to  demodulate  the  measured 
spectrum  to  remove  the  effect  of  the  surface  reflection  In  the  frequency 
domain  have  been  made'*  These  efforts  have  generally  been  hampered  by 
the  lack  of  a  realistic  model  for  the  surface  reflected  pressure.  The 
reflected  pressure  history  is  not  always,  as  is  often  assumed,  simply  the 
mirror  image  of  the  direct  pressures  with  an  appropriate  time  delay.  As 
pointed  out  by  Parkins',  scattering  at  the  interface  tends  to  destroy  the 
coherency  of  surface  reflected  arrivals.  The  noise  introduced  by  the 
scattering  seriously  erodes  the  simple  ideal  modulation  when  surface  rough¬ 
ness  Is  important.  A  more  accurate  model  of  the  surface  reflection  process, 
taking  surface  scattering  into  account,  would  be  needed  for  this  method  to 
prove  more  generally  applicable.  In  some  instances,  however,  the  simple  model 
may  be  acceptable  (for  example,  see  paragraph  7.6  and  reference  15). 


♦  The  300- and  500- ft.  curves  of  reference  2  are  not  based  entirely  on  actual  data 
at  those  depths.  The  curves  are  Interpolated  between  the  deeper  curves  and 
Stockhausen's  shock  were  data  (reference  10 )  at  high  frequencies  and  extrapolated 
from  the  deeper  curves  at  low  frequencies. 
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4.3*5  Efforts  to  eliminate  the  surface  reflection  from  measured 
pressure  records  by  digitally  processing  them  and  setting  all  negative  pressures 
equal  to  zero,  as  done  in  reference  13  Is  unrealistic  in  that  the  energy 
contribution  of  the  negative  pressure  phases  nt  direct  arrival,  vhich  gay 
be  appreciable  et  lov  frequencies,  is  eliminated. 

4.4  Of  the  above  methods,  the  use  of  idealized  free  water  pressure-time 
curves  appears  to  be  the  most  promising  for  the  determination  of  shallow 
explosive  source  levels.  In  the  method  of  the  present  report,  the  best 
available  Information  on  explosive  pressure-time  histories  has  been  utilized  to 
construct  qua s i -the or et i cal  pressure-time  curves  for  several  shallow  explosion 
geometries  of  interest.  We  feel  that  these  idealized  pressure  histories 
represent  an  Improvement  over  those  methods  previously  used.  The  details  of 
their  construction  are  given  in  the  following  section. 

r>.  CONSTRUCTION  OF  IDEALIZED  FRESSUUt-TIME  HISTORIES 


5.1  Idealized  pressure-time  histories  for  severed  combinations  of 
charge  weight  and  depth  have  been  constructed  utilizing  the  best  available 
information  for  deeper  explosions  and  extrapolating  to  the  shallower  depths. 
The  details  of  this  procedure  are  given  below. 

The  major  body  of  information  utilized  was  reported  by  Slifko^. 
This  report  documents  the  analysis  of  a  large  body  of  data  gathered  by  the 
Naval  Ordnance  Laboratory  (NOL).  A  total  of  5 6  charges  weighing  1,  8  and 
50  lb  were  fired  et  depths  from  500  to  14,000  feet.  Pressure-time  records 
were  obtained  with  an  LC32  hydrophone  located  directly  above  the  charge  and 
185  feet  from  the  surface.  Seml-emplrlcal  equations  were  derived  vhlch  give 
the  various  parameters  of  the  pressure-time  record  as  functions  of  depth, 
range  and  charge  weight.  Theoretical  considerations  were  utilized  to 
separate  the  effects  of  range  and  charge  depth. 

5.2  In  addition,  information  from  other  available  sources  was 
Incorporated  Into  the  theoretical  curves  used  for  this  study.  The  various 
parameters  used  as  inputs  to  the  construction  of  the  idealized  pressure 
histories  are  indicated  in  Figure  2.  Table  2  gives  the  functional  forms 
and  sources  far  these  parameters. 
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Table  2 


FUNCTIONAL  DEPENDENCE  OF  PRESSURE-TIME  PARAMETERS 


•aaeter 

Functional  Dependence 

Source 

Ref 

P 

O 

2.08  x  10*  (w1^)1'13 

Slifko 

16 

Pmin 

-77  Z1//,3Wly/3/R 

Slifko 

16 

F1 

3300  Wly/3/R 

Slifko 

16 

P2 

.  22  Fx 

Unpublished 
data  at  NOL 

P3 

.10  P 

Unpublished 
data  at  NOL 

P4 

•  03  PL 

Unpublished 
data  at  NOL 

9 

5.8  x  10" V^3  (VL/3/R) 

-0.22 

Arons 

18 

T 

PP 

T 

np 

r 

•  555 

Slifko 

16 

-  KW1^3.  K  • 

176 

z 

3.10 

Slifko . 

16 

Tbp 

1.10 

Slifko 

16 

T1 

► 

< 

4.34 

Slifko 

16 

T2 

3.06 

Arons, 

Slifko,  and 
Carter 

17 

T3 

2.48 

Arons, 

Slifko,  and 
Carter 

17 

T4  . 

\ 

2.31 

Arons, 

Slifko,  and 
Carter 

17 

All  pressures  In  psl  and  times  are  in  seconds 

R  -  slnnt  range,  ft. 
w  ’  charge  weight,  lbs. 

-  *  charge  depth  +  33#  ft* 
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*>•3  The  shock  wave  was  constructed  lnttJall^  exponential  with  a  decay  rate,  9, 
as  given  hy  the  similitude  equation  of  Arons^".  The  exponential  decay  was  continued 
out  to  approximately  29.  The  impulse  In  the  first  positive  phase  (shock  wave) 
anu  in  the  first  bubble  pulse  positive  phase  was  calculated  from  the  relations 
givsr.  by  Giifku--.  We  had  no  information  on  tne  impulse  in  the  later  bubble 
pulses  and  the  minimum  pressure  of  the  later  negative  phases.  From  Inspection 
of  unjaibllshed  data.  It  was  found  that  the  maximum  underpressure  in  the  successive 
negative  phases  tended  to  become  smaller.  These  parameters  did  not  completely 
define  the  pressure  signature,  but  it  was  felt  that  the  uncertainties  in  the  later 
states  would  not  be  of  primary  importance.  The  pressure-time  curves  were  hand 
drawn  according  to  the  information  given  in  Table  2  as  described  above.  The 
records  were  digitized  by  use  of  the  Telereadex  and  were  analyzed  digitally  using 
NOL'e  CDC  6100  computer. 

o.l  The  first  stage  of  the  digital  analysis  was  checking  the  digitizea 
waveform  for  compliance  with  the  information  inputs.  Minor  modifications  were 
made  to  the  first  and  second  positive  phases  to  bring  the  Impulse  in  these  phases 
to  within  1056  of  those  calculated  from  Slifko's  relationships.  The  negative 
phases  were  modified  so  as  to  produce  near  zero  residual  impulse  out  to  the  fourth 
bubble  pulse. 

o.5  The  analysis  was  performed  using  a  system  of  computer  codes  developed 
at  tfOL  and  known  by  the  acronym,  MR.  W1SARD.  The  package  is  described  in  detail 
In  Reference  19>  Using  MR.  WISARD,  the  spectra  were  computed  using  the  fast 
Fourier  transform  technique  of  Cooley  and  Tukey20,  with  a  bandwidth  of  about 
1  Hz.  The  spectra  were  also  integrated  in  octave  and  l/3 -octave  bands  (the  bands 
used  are  indicated  in  Table  3)  to  facilitate  comparison  with  spectra  determined 
by  analogue  methods.  The  filter  bands  used  were  ideal  (not  physically  realizable} 
in  that  they  passed  all  the  energy  between  the  limiting  frequencies  and  no  energy 
from  outside  those  frequencies. 

5.6  In  order  to  make  comparisons  with  measured  data  for  the  shallow  conditions 
an  idealized  surface  reflection  was  introduced  for  some  parts  of  the  analysis. 

The  reflection  was  assumed  to  cause  a  180°  phase  shift  with  no  losses  at  the 
interface.  The  amplitude  of  the  reflected  pulse  was  adjusted  for  the  difference^ 
in  travel  distance  between  it  and  the  direct  arrival.  As  suggested  by  Christ isu  -  , 
the  surface  reflection  was  put  in  as  a  modification  in  the  time  domain,  rather 
than  a  modulation  of  the  energy  spectrum.  In  our  particular  computational  scheme, 
the  effort  required  to  modify  the  tine  domain  functions  was  far  less  than  that 
needed  to  correct  for  surface  reflections  In  the  frequency  domain. 

i>.7  These  Idealized  pressure- dice  records  represent  the  best  available 
synthesis  of  empirical  sund  theoretical  knowledge  available  to  the  authors  at  this 
time.  The  extrapolation  from  deeper  to  relatively  shallow  depths  represents  a 
major  uncertainty.  In  the  absence  of  a  good  theoretical  treatment  of  shallow 
explosions,  it  was  felt  that  the  above  represented  the  most  promising  approach 
to  the  problem.  The  results  and  comparisons  of  the  following  two  sections  shed 
light  on  the  validity  of  the  entire  procedure. 
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6.  RES  UUTS 

6.1  Idealized  pressure-time  curves  have  been  constructed  using  the  methods 
detailed  in  the  previous  section  for  several  conditions  of  interest.  The  MK  6l 
signal  np.d-rv-tcr  Sound  (S'JE)  charge  contain*  l.S  lb  of  a^luelve  natefial  ("SI) 
and  is  detonable  at  60  ft.  and  8 00  ft.  The  SUB  MK  57-0  and  MK  82-0  also  contain 
1.8  lb  of  TlfT  and  have  nominal  firing  depths  of  8 00  ft.  and  60  or  300  ft.,  respec¬ 
tively.  To  obtain  source  energy  levels  for  these  devices,  pressure-tine  records 
for  1.6  lb  of  TlfT  detonated  at  depths  of  60,  300  and  800  ft.  vere  constructed. 

For  the  60-foot  condition,  no  direct  measurement  of  the  source  level  is  possible. 

For  the  300-foot,  1.8-lb  source,  with  a  100-yard  reference  range,  the  angvii 
position  of  the  receiver  may  vary  the  arrival  time  of  the  surface  reflection  from 
being  coincident  with  the  direct  arrival  directly  above  the  charge  at  tha  reference 
range  t  >  being  between  the  third  and  fourth  bubble  pulses  directly  below  the 
charge.  Hence,  direct  measurement  of  the  source  levels  will  generally  appear  to 
exhibit  directivity.  The  deeper  800-foot,  1.8-lb  condition  offers  no  problem 
since  the  direct  and  surface  reflected  arrivals  are  easily  time  separable.  This 
condition  Is  Included  for  completeness  and  for  comparison  purposes. 

6.2  In  addition  to  the  1.8-lb  charge  weight,  a  weight  of  3  lb  was  also 
treated  at  depths  of  detonation  of  60  and  500  feet.  These  conditions  vere  Included 
because  they  have  been  used  in  important  oceanic  acoustic  experiments.  The  500-foot, 
3  lb  source  vfoich  offers  no  apparent  source  directivity  is  included  for  completeness, 
while  the  60-foot,  3-lb  source  oust  be  treated  by  the  method  developed  here  for 
shallow  charges. 

6. 3  To  summarise,  the  five  conditions  calculated  were  1.8  lb  at  60,  300,  and 
800- foot  depths  end  3  lbs  at  60  and  500- foot  depths.  Two  representative  pressure- 
time  curves  for  3  lb  at  60  feet  and  1.8  lb  at  800  fset,  are  shown  in  Figures  3 
and  u. 


6.14  Energy  spectra  were  computed  for  the  five  conditions  and  are  shown  in 
Figures  5-9.  These  spectra  were  computed  with  a  resolution  bandwidth  of 
approximately  1  Hz.  Each  spectrum  shows  a  complex  pattern  of  peaks  and  nulls,  the 
primary  features  of  lAlch  may  be  related  easily  to  the  time  domain  characteristics. 
The  highest,  peak  occurs  at  the  bubble  fundamental  frequency*,  and  succeeding  peaks 
at  approximately  higher  harmonics  of  this  frequency.  The  dominance  of  the  first 
bubble  pulse  in  the  shape  of  the  spectrum  is  due  to  the  fact  that  tha  first  bubble 
pulse  contains  significantly  more  energy  than  the  following  pulses.  The  later 
pulses  do  have  am  effect  on  the  spectrum,  as  they  introduce  secondary  Interference 
patterns  (tdilch  are  conspicuous  at  the  lower  harmonics)  and  shift  the  peak  locations 
at  higher  frequencies. 

6. 5  Although  the  source  level  at  a  range  of  100  yards  is  fully  defined  by  the 
detailed  oar row-band  spectra  of  Figures  5-9,  many  workers  in  ocean  acoustics  use 
frequency  bands  of  l/3-octave  and  1  octave  in  their  analyses.  We  have,  therefore, 
Integrated  the  energy  spectra  in  these  energy  bands,  and  present  the  results  below. 

6.6  Table  3  gives  the  energy  distribution  for  l/3-octave  bands  and  Table  b 
gives  it  for  the  octave  bands  for  the  five  conditions  that  have  been  studied.  These 

*  The  reciprocal  of  the  time  of  the  first  bubble  pulse. 
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vmluss  are  graphically  displayed  in  Figures  10-11.  1%e  energy  rises  from  low 

values  at  low  frequency,  to  a  maxi bub  near  the  bubble  pulse  fundaaentel  frequency. 

It  falls  off  thereafter,  and  at  higher  frequencies  approaches  the  6  db/octave 
fell off  predicted  for  exponential  waves.  This  is  due  to  the  dominance  of  the 
shock  an  at  high  frequencies.  The  relative  smoothness  of  the  octave  hand  energy 
is  due  +.n  the  fact  that  the  intogiibiua  bands  are  generally  vide  enough  to  average 
out  the  peaks  and  nulls  present  in  the  detailed  spectrum  (l  Hi  bandwidth* )  of 
Figures  5-9.  The  l/3  octave  analysis  shows  soae  Indication  of  the  detailed 
structure  of  the  spectrum.  It  should  be  pointed  out  that  the  apparent  irregularity 
of  the  l/3 -octave  spectra  does  not  represent  spurious  behavior.  The  lav  frequency 
Analysis  bands  are  generally  narrow  with  respect  to  the  features  of  the  detailed 
spectra,  and  nay  contain  a  paak,  a  null  or  any  coablnatlon  of  peaks  and  nulls. 
Therefore,  the  energy  nay  change  in  a  seemingly  sporadic  way  from  band  to  band. 

At  higher  frequencies,  as  the  bandB  become  wider,  they  encompass  a  number  of  peakc 
end  nulls  of  the  spectrum.  When  the  energy  is  Integrated  in  these  bands,  the 
abrupt  chsuigeo  are  smoothed  out,  and  the  spectrum  approaches  the  smooth  behavior 
of  the  octave  spectra.  Hence,  the  apparent  Irregularity  of  the  spectra  in 
l/3-octava  bands  is  seen  to  be  caused  by  the  relationship  of  the  analysis  bandt 
uted  to  the  features  of  the  detailed  spectrum.  The  fluctuations  are  real,  ana 
represent  genuine  fluctuation  in  the  source  level  with  frequency. 

6.7  In  the  figures  showing  the  l/3-and  octave  band  spectra,  the  data  points 
have  been  connected  with  straight  lines.  This  has  been  done  for  convenience  in 
demonstreting  the  variation  in  energy  from  band  to  band.  It  is  import  suit  to  note 
that  these  figures  may  not  generally  be  used  to  interpolate  between  bands.  For 
example,  suppose  we  wish  to  determine  the  source  level  for  a  60-* t  3-lb  TlfT  charge 
at  25  Hz.  The  detailed  spectrum  of  Figure  3  shows  a  null  at  this  frequency,  with 
energy  of  about  b  db.  The  l/3-octave  energy  calculation  at  25  Hz  in  Figure  10  Is 
19.5  db.  Both  these  numbers  are  correct,  the  former  giving  the  energy  mt  25  Hz 
for  an  approximate  1  Hz  bandwidth,  and  the  latter,  the  average  energy  in  a 

l/3 -octave  band  with  a  geometric  mean  frequency  of  25  Hz.  These  numbers  would  be 
expected  to  be  different,  due  to  the  difference  In  the  physical  quantities 
represented.  Again,  though  they  differ  greatly,  both  numbers  are  valid.  In  regions 
of  relatively  slov  variation  in  the  octave  spectrum,  say  for  frequencies  at  least 
twice  the  bubble  fundamental,  linear  interpolation  may  give  an  indication  of  the 
energy  to  be  expected  in  an  octave  band  centered  at  the  glvan  frequency.  This 
would  only  be  a  rough  approximation.  Linear  interpolation  at  lower  frequencies, 
however,  would  be  extremely  dubious. 

6.8  If  one  were  using  the  octave  spectrum  of  Figure  10  In  an  attempt  to 
determine  the  energy  at  25  Hz — regardless  of  Aether  he  is  interested  in  a  band  of 
1-Rz,  1-octave  or  l/3-octave,  he  would  not  be  able  to  extract  from  this  figure  the 
information  desired.  He  slight  erroneously  Interpolate  linearly  between  the  bands 
centered  at  16  Hz  and  31. 3  Hz  to  obtain  21  db— a  value  having  no  real  meaning. 

From  the  exercise  in  this  and  the  previous  paragraph,  three  different  values  for 
energy  et  25  Hz  have  been  obtained.  The  first  two,  U  db  and  19.5  db,  refer  to  the 
energy  in  a  1-Hz  band  centered  at  25  Hz  and  the  energy  in  a  l/3-octave  hand 
centered  at  25  Hz,  respectively.  The  21  db  at  23  Hz  of  this  paragraph  has  no 
correct  interpretation.  It  Is  clear  that  the  source  level  for  a  particular 
application  must  be  carefully  determined  and  interpreted  in  terms  of  the  appropriate 
bandwidth  analysis  used. 
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TABLE  '■ 

SOURCE  LEVELS  IK  OCTAVE  BANDS-ENERGY  IK  db  re  I  erg/cn2/Hz 
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6.9  Fffect  of  Small  Variation*  in  Burst  Depth 

6.9.1  It  is  a  fact  of  experimental  life  that  explosions  often  do  not 
detonate  at  precisely  their  nominal  depths.  The  depth  span  to  be  expected  for  SUS 
charges  are  £0-j0  feet  for  a  nominal  oO-ft  bur at ^1.  end  700-900  feet  ?at  ■ 

300- ft  Burst^'^J .  Approximately  the  ssuae  percentage  variation  In  buret  depth 
would  probably  occur  with  other  nominal  depth  settings.  These  seemingly  small 
variations  in  burst  depth  may  have  a  significant  effect  on  the  source  level  in 
certain  of  the  l/3*octave  and  octave  bands  in  the  tine  domain,  the  primary  effect 
of  depth  variations  is  a  change  in  the  bubble  periods  and  other  characteristic 
tines,  vhich  all  vary  as  Z"5/°,  uhere  Z  is  the  charge  depth  +  33  ft,  (see  Table  2). 
The  variation  of  the  minlnum  pressure  in  the  first  negative  phase  Is  considered  a 
second  order  effect. 

6.9.?  In  the  frequency  domain,  the  location  of  the  peak  at  the  bubble 
fundamental  is  shifted,  and  the  energy  at  the  peak  decreases  with  increasing  burst 
depth,  figure  15  shovs  the  detailed  spectra  far  1.8  lbs  of  TOT  at  burst  depths  of 
700,  800,  and  900  feet.  The  bubble  fundamental  frequencies  for  these  depths  are 
b-6,  51,  and  56  Hz,  respectively.  The  pattern  of  variation  around  the  first  peak 
in  the  spectrum  is  generally  repeated  at  the  succeeding  peaks,  but  the  variation  is 
increasingly  masked  by  the  secondary  Interference  patterns.  The  consistent 
variation  of  energy  with  burst  depth  near  the  bubble  fundamental  may  be  discerned 
in  the  l/3-octave  band  analysis.  The  differences  In  energy,  AE,  between  the  source 
level  at  8 00  feet  and  that  700,  750,  850,  and  900  feet  are  tabulated  In  Table  5- 
Up  to  about  50  Hz,  the  l/3-octave  band  levels  display  a  fairly  consistent  variation 
with  depth,  the  energy  level  decreasing  for  increasing  burst  depth.  Above  this 
frequency,  the  pattern  of  variation  seems  store  random.  Table  6  gives  a  similar 
tabulation  for  a  nominal  60-foot  burst  of  a  1.8-lb  TOT  charge.  The  bubble 
fundamental  frequency  (-£  Hz  for  a  60-ft  depth)  Is  too  low  to  be  of  practical 
importance  in  most  applications.  The  variations  depicted,  then,  are  all  for 
frequencies  considerably  above  the  fundamental,  and  have  the  same  apparently  random 
character  as  the  higher  frequency  variations  of  Table  5. 

6.9.3  The  variation  In  source  level  caused  by  the  above  variations  In  depth 
may  be  as  high  as  6.1  db  (in  the  20  Hz  l/3-octave  band)  for  the  800  foot  nominal 
depth,  and  as  much  as  2.5  db  (in  the  31* 5  Hz  l/3  octave  band )  for  the  60-foot 
nominal  depth.  These  variations  may  be  quite  significant  In  some  applications. 

The  Implication  here  is  that  the  depth  of  burst  must  be  accurately  determined  in 
order  for  the  source  level  to  be  given  with  reasonable  accuracy.  Merely  using  the 
nominal  depth  may  lead  to  considerable,  seemingly  disproportionate,  errors 
(6  db  represents  a  factor  of  h  in  energy). 

6. 10  Surface  Reflections 

6.10.1  As  mentioned  above,  apparently  some  workers  have  erroneously  Included 
surface  reflections  as  part  of  the  source  level.  The  effect  of  apface  reflections 
on  the  energy  spectrum  has  been  described  in  detail  by  Christian  .  Briefly,  the 
effect  is  to  impose  a  modulation,  due  to  the  reflection,  on  the  spectrum  of  the 
direct  arrival. 

6.10.2  As  an  example  of  this,  take  the  case  of  a  3-lb  TOT  charge  detonated  at 
50  feet.  The  maximum  time  delay  between  the  direct  and  surface  reflected  arrivals 
is  2h  msec.  This  maximum  delay  Is  obtained  when  the  measurement  13  made  directly 
beneath  the  charge.  Since  this  maxinum  reflection  time  is  shorter  than  the  length 
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TABLE  : 


DEPTH  VARIATION  FOP  AN  8 OO-FOOT,  1,8-LB  SOURCE 

AE  (E.  -  eqqq)  (<3b  re  1  erg/cin'/Hz) 

_ _ Range  100  yds 

l/3  Octave  Band  Source  Depth 

Center  Frequency _ 700*  .  75b1  0501  900 1 


12.5 

1.8 

.9 

-.7 

-1.4 

16 

2.8 

1.3 

-1.1 

-2.0 

20 

3-3 

1.5 

-1.5 

-2.3 

25 

2.8 

1.4 

-1.4 

-2.7 

31.5 

2.2 

1.1 

-1.0 

„p.n 

40 

2.4 

1.2 

-1.0 

-1.3 

50 

1.2 

r» 

•  1 

-.9 

-1.7 

63 

1.2 

-.4 

-.1 

O.C 

80 

2.0 

-1.0 

.1 

o.1. 

100 

0.0 

.1 

-.2 

-  •  5 

125 

1.7 

-1.3 

.8 

1.1 

160 

2.2 

-.4 

-.2 

-1.2 

200 

1.0 

-.5 

■1.6 

-2.3 

250 

1.2 

.1 

0.0 

•  3 

315 

1.6 

1.7 

.8 

1.1 

400 

.5 

.4 

-.4 

-.3 

500 

.4 

.3 

-.1 

-.5 

630 

.  6 

.4 

-.1 

-.2 

800 

.8 

.5 

.2 

-.1 

1000 

.7 

.5 

.2 

0.0 
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TABLE  6 


nrOTU  If  A  D  7  iTTAk'  .  - _ _ 

— . A  uak  A  wruui,  J.«0-laB  £>UUftU£ 


AS  (Ed  -  n^)  {db  re  1  •rg/cm2/H*) 


Range  100  yds 


l/3-0etave  Bend 
Center  Frequency 

12.5 
16 
20 
25 

31.5 
1+0 
50 

63 

80 

100 

125 

160 

200 

250 

315 

400 

500 

630 

800 

1000 


Source  Depth 

W  55 1  r  fef 


.4 

.4 

•5 

.0 

-.7 

.2 

.0 

-.6 

1.2 

-1.2 

.1 

-.4 

1.1 

1.5 

-.5 

-1.0 

1.1 

.6 

w 

.7 

.4 

-.1 

-1.8 

-1.1 

.0 

1.5 

-1.7 

.3 

-.2 

.2 

.0 

-•3 

*•3 

-.4 

*.3 

-.2 

-.1 

.1 

-•3 
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of  the  pulse  train,  the  reflection  is  not  separable  from  the  direct  arrival. 

Figure  16  shows  the  detailed  spectra  for  the  ideal  pulse  with  no  reflection  and 
with  a  reflection  at  24  msec.  As  nay  be  seen,  the  reflection  imposes  a  modulation 
on  the  free  water  spectrua,  vith  nulls  at  approximately  integral  multiples  of  the 
inverse  of  the  surface  rejection  time  ,  1  >.  -  „  ihe  effect  of  this 

l24  msec  "  41,7 

modulation  on  the  l/3- octave  band  spectra  is  shown  in  Figure  17.  The  errors 
introduced  by  including  the  surface  reflection  as  part  of  the  source  level  are 
large  and  seemingly  random.  For  example,  the  surface  reflection  adds  about  4  db  at 
20  Hz,  and  subtracts  about  5  db  at  40  Hz.  lfo  simple  correction  of  a  constant 
number  of  db  will  suffice.  This  example  Indicates  the  primary  importance  of 
eliminating  the  surface  reflection  vfcen  a  source  level  la  to  he  obtained  from  a 
treasured  pressure-time  history. 

6.11  In  the  following  section,  some  preliminary  comparisons  of  the  results 
of  the  present  study  with  previous  calculations  and  measured  data  are  presented. 

7.  COMPARISONS  WITH  PREVIOUS  WORK 


7.1  In  order  to  establish  the  validity  of  the  shallow  (60  foot)  source  levels 
of  the  present  study,  the  theoretical  model  should  he  compared  vith  appropriate 
measurements.  Of  course,  a  surface  reflection  would  have  to  be  added  to  the 
idealized  pulse  in  order  to  make  such  comparisons.  For  such  a  comparison  to  be 
realistic,  however,  experimental  conditions  (e.g.,  accurate  charge  and  receiver 
locations,  equipment  response  characteristics)  must  be  known  in  greater  detail  than 
is  reported  for  available  published  data.  Efforts  to  obtain  existing  unpublished 
data  for  analysis  and  comparison  have  been  unsuccessful  up  to  the  time  of  this 
writing. 


7.2  In  the  absence  of  such  shallow  data,  comparisons  with  deeper  date  have 
bean  made.  These  comparisons  verify  the  accuracy  of  the  method  used  in  this  study, 
except  (as  indicated  above )  for  the  vital  extrapolation  of  seml-emplrlcal  pressure 
time  relations  from  deep  (>500  feet)  to  shallow  depths,  as  discussed  in  section  5. 
The  extrapolation  has  been  partially  evaluated  by  comparisons  with  data  at  a  300 
foot  charge  depth,  but  further  comparisons  will  be  necessary. 

2 

7.3  The  present  results  are  compared  with  those  of  Christian  in  Figures 
18  through  20.  The  3°0  and  5 00  foot  curves  are  from  reference  2,  ifcile  the  800 
foot  source  level  curve  is  previously  unpublished  data25.  The  idealized  source 
spectra  have  been  integrated  in  the  same  octave  beads  and  scaled  in  the  same 
manner  as  done  by  Christian.  The  maximum  difference  shown  on  these  three  figures 

is  —  2.5  <2b  with  a  mean  absolute  difference  of  1  db.  Considering  that  Christian's 
300-and  500-foot  curves  are  extrapolations  that  apparently  fall  to  allow 
sufficient  energy  contributions  from  the  bubble  phases*,  the  results  dereived  from 
this  study  are  believed  to  be  more  fcorrect  for  these  shallower  conditions  than  the 
published  ones.  Agreement  between  the  present  study  and  Christian's  800-foot  data 
is  excellent  (Figure  20 ). 

12 

7.4  The  data  of  Turner  and  Scrimger  offer  further  valuable  comparisons. 

Their  source  level  data  for  1  lb  of  pentolite  at  depths  of  295  and  555  feet  have 
been  compared  vith  the  Idealized  results  of  this  study,  and  this  comparison  is 


*  See  footnote  on  page  6. 
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presented  in  Table  7.  The  energy  spectra  of  the  idealized  pressure  records  were 
computed  after  adjusting  the  pressure  and  time  scales  to  those  of  a  1-lb  TNT  charge 
at  the  same  ranges  and  depths.  No  allowance  for  the  different  explosive 
compositions  was  included  in  the  calculation.  For  both  test  conditions.  Turner 
and  Scri«jr«*r  «hmj  energy  nt  lov  fro^ucuciee  and  less  at  nigh  frequencies  than 

the  present  results.  This  is  the  same  general  pattern  as  shown  In  their  published 
comparison  with  Christian's  source  level  curves,  though  the  magnitude  of  the 
differences  cure  reduced  in  the  present  study.  At  295  feet,  the  aaxlmum  deviation 
in  Table  7  is  2.6  db.  (The  mean  absolute  deviation  is  1.6  db. )  This  is 
considerably  better  than  the  555-foot  comparison,  ldiich  has  a  maximum  deviation  of 
nearly  .?  db.  (The  mean  absolute  deviation  is  2.6  lb. )  Part  of  these  differences 
may  be  attributable  to  the  fact  that  pentollte  has  about  20$  more  shock  wave 
energy  than  TNT  but  the  major  differences,  especially  at  555  feet,  are  unexplained. 
Differences  in  processing  techniques,  euch  as  recording,  range  correction, 
digitization  and  computer  analysis  techniques  are  possible  sources  of  these  large 
deviations. 

7.5  Several  additional  samples  of  data  available  at  NOL  have  been  analyzed 
for  comparison  with  the  idealized  pectre.  These  data 'are  part  of  another  ■ 
collection  reported  by  Christian  and  others  in  a  classified  report  and  will  be 
referred  to  as  Squaw  3hots.  Three  samples  of  data  from  charges  at  a  nominal  depth 
of  900  feet  are  compared  with  the  present  results  in  Figure  21.  The  actual  depths 
were  835,  967,  and  676  feet.  The  recordings  were  made  at  ranges  of  1200,  1000, 
and  266  yards  respectively  from  the  three  shot,?.  The  sources  were  Mk  6l  SUS  charges 
(l.3  lb  of  TNT).  The  three  records  were  scaled  ho  an  800- foot  source  depth  by 
adjusting  the  time  3cale;  spherical  spreading  vas  used  to  correct  to  a  100-yard 
reference  range.  The  energy  was  computed  in  }./ 3-octave  bands.  The  source  levels 
for  these  three  Squaw  shots  and  for  the  idealized  800-foot  curve  (in  l/3-octave 
bands)  are  shown  in  Figure  21.  Note  that  the  scatter  among  the  three  Squaw  shotr 
is  on  the  order  of  1  to  2  db.  The  differences  between  the  Squaw  shots  and  the 
idealized  curve  sire  generally  less  than  2  db  up  to  about  5°0  Hz.  The  differences 
above  500  Hz  have  not  been  resolved  at  this  time,  but  cmr  assumption  of  simple 
spherical  spreading  to  correct  for  range  differences  probably  accounts  for  them. 

The  general  pattern  of  irregularity  shown  in  the  ii*<.ta  is  well  modelled  by  the 
idealized  curve.  The  dashed  curve  also  shown  in  Figuze21  is  discussed  in  section 
7.9  below. 

7-6  A  fourth  Squaw  shot  was  selected  to  provide  a  mor*  shallow  comparison. 

This  was  another  Mk  61  SUR  detonated  at  300  feet.  The  broau  her  I  o"*ssure-time 
record  shows  a  surface  reflected  arrival  at  2b. 7  msec  after  .  -ect  arrival. 

For  comparison,  a  180*  phase  shifted  surface  reflection  was  iai.  u*.vt)d  into  the 
idealized  pressure-time  curve  at  2^.7  msec.  Both  theoretical  and  measured  pressure 
histories  were  analyzed  in  l/3-octave  bands  and  are  shown  in  Figure  22.  The 
agreement  again  is  good;  the  general  pattern  of  peaks  and  nulls  is  especially  close. 
Considering  the  crude  nature  of  the  surface  reflection  model  used,  the  agreement  is 
excellent.  This  is  an  indication  that  the  extrapolation  of  our  idealized  waveform 
parameters  from  >  50C  to  300- ft  source  depth  is  a  valid  procedure.  'Hie  further 
extrapolation,  to  a  60-foot  source  depth,  remains  uncertain. 

7.7  Since  the  calculation  of  Weston  has  long  been  the  standard  one  for  source 
levels,  a  comparison  between  the  present  results  and  those  of  Weston  is  included 
here.  Before  presenting  this  comparison,  some  comments  regarding  the  differences 
between  the  two  studies  are  in  order. 
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TABLE  7 

COMPARISON  OF  TURNER  AND  SCRDttHt’S  RESULTS  (Ref.  12)  WITH  THE 


(Energy  In  db  re  1  erg/cm^/HE) 

Range  100  yds 

1  lb  at  295*  1  lb  at  555* 


Frea.  Bend 

Ref  12 

Present  Study 

A 

5SU£ 

Present  Study 

A 

20-40 

17.6 

16.4 

1.2 

13.0 

8.7 

4.3 

40-80 

14.7 

12.5 

2.2 

17.7 

14.0 

3.7 

80-160 

12.4 

10.1 

2.3 

13.6 

9.0 

4.6 

160-370 

9  A 

6.8 

2.6 

9.9, 

6.2 

3.7 

370-640 

6.0 

4.1 

1.9 

6.4 

4.3 

2.1 

640-1280 

-.9 

•  3 

-1.2 

-.5 

1.2 

-1.7 

1280-2560 

-4.8 

-3.6 

-1.2 

-4.8 

-3.4 

-1.4 

2560-5120 

-10.1 

-8.7 

-1.4 

-10.0 

-8.6 

-1.4 

5120-10240 

-15.0 

-14.4 

-.6 

-14.9 

-14.4 

-.5 

Mean  absolute 

deviation 

1.6 

2.6 
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7.8  As  mentioned  Above,  in  the  low  frequency  region,  Weston  uses  an  Impulse 
formula  to  calculate  the  energy  spectrum.  Both  the  shock  vave  positive  Impulse 
and  that  of  the  first  two  bubble  pulses  are  included,  with  a  negative  impulse  equal 
to  their  suu  added  to  achieve  M'c  residual  total  impulse.  The  formulas  used  far 
the  Impulse  vers  the  bssi  available  at  that  time,  but  must  be  examined  carefully 
in  the  light  of  subsequent  work.  Weston  uses  the  shock  wave  positive  Impulse 
formula  reported  by  Arons^®.  This  formula  was  synthesized  from  experimental 
■•tjcords  of  large  harges  (2>-76  lb)  fired  in  shallow  water.  Since  the  time  of 
arrival  of  tha  surface  reflection  was  very  short  in  relation  to  the  duration  of 
the  shock  wave  positive  phase,  the  impulse  could  not  be  determined  for  the  entire 
shock  wave  phase  of  the  pulse.  Although  the  point  is  not  made  clear  by  Arons, 
comparison  of  the  fonnula  with  that  reported  by  Aro.-.c  in  reference  26  indicates 
that  the  impulse  represented  by  Arons'  formula  is  integrated  out  to  only  6. 7ft ^  as 
was  •uctooary  at  that  time.  Arms  and  .eunie2?  indicate  that  the  impulse  at  6. 7P 
i 3  some  hing  less  than  63$  of  the  totil  impulse  of  the  shock  wave  for  a  250-foot 
deep  chit.  In  addition,  the  results  of  Arons  show  no  depth  variation  in  the 
impulse,  since  all  of  Arons'  data  were  from  approximately  the  same  depth.  The 
shock  wave  impulse  formula  of  Slifko,  as  used  In  this  study,  represents  the  total 
impulse  of  the  first  positive  phase  (out  to  the  first  zero  pressure  crossing),  and 
indicates  a  significant  depth  variation.  The  ratio  of  Slifko's  to  Arons'  impulse 
formula  as  a  function  of  charge  depth  is  shown  in  Figure  23**  Shallower  than  500 
feet,  Slifko's  formula  represents  an  extrapolation  from  his  measured  data.  Ibis 
extrapolation,  while  uncertain,  is  probably  a  more  reasonable  indication  of  the 
Impulse  at  shallow  depths  than  Arons'  formula  on  two  counts: 

1.  Arons'  formula  represents  the  impulse  out  to  only  6.7°  While  flifko's 
is  for  the  entire  shock  wave  p>aitive  phase,  and 

2.  Sllfko  takas  depth  variation  into  account  while  Arons  does  not. 

Aron's  and  Slifko's  famtlas  agree  mt  >00- ft  depth.  At  more  shallow  depths, 

Slifko  gives  more  impulse;  at  deeper  deptha,  less  impulea.  The  impulse 
formula  used  by  Weston  for  the  first  bubble  pulsa  Impulse  is  also  quite 
different  fro*  that  used  in  the  present  study**.  One  would  expect 
considerable  differences  between  Weston's  analysis  and  the  present  study 

In  tha  low  frequency  region.  It  is  felt  that  tha  waveforms  of  the  present 
study  allow  a  more  complete  and  accurate  representation  of  spectral  energy 
content  than  the  approximations  of  Weston. 

7*9  In  order  to  make  comparisons,  Weston's  analytical  calculation  has  been 
repeated  for  the  1.8-lb  SUE  charge  at  depths  of  60  and  800  feet.  The  results  have 
been  Integrated  in  l/3-octave  bands  and  are  compared  with  the  corresponding  present 
results.  The  1/3-octave  band  differences  are  indicated  in  Table  8;  the  detailed 
spectra  are  compered  in  Figures  24  and  25 >  For  1.8-lb  at  60  ft  condition,  the 
present  results  indicate  considerably  more  energy--as  much  as  8  db  in  the  very  low 
frequency  bands.  This  is  largely  due  to  the  different  impulse  formulas  used.  In 
high  frequency  bands,  Weston  shows  2  to  3  db  more  energy  for  the  60-ft  source. 

For  an  800- ft  source  Weston's  formulation  gives  much  more  energy  at  very  low 
frequencies— again  due  largely  to  the  different  impulse  formulas;  at  higher 
frequencies,  Weston  shows  1  to  3  db  more  energy.  The  differences  for  both  burst 
depths  are  significant.  To  resolve  this  discrepancy,  comparisons  with  data  are 
needed.  As  noted  above,  we  have  no  suitable  60- ft.  data  for  comparisons.  For  the 

-  703 

*  The  factor  (-tk)  has  been  neglected  In  this  ratio. 

**  Wegton's  equation  7  contains  a  typographical  error.  The  exponent  of  W  should  be 
2/7. 
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TABLE  8 


l/3-OCTAVE  BAND  COMPARISON  OP  CALCUIATIOIB  USING  WESTON'S 

ANALYTICAL  FORMULATION  WITH  THOSE  OF  THE  PRESENT  STUDY 

Range  100  yds  „ 

AE  (] Resent  reaulte  '  Keeton  ln  db  ™  1  er«/cM  **) 
l/3-Oct eve  band 

center  frequency  1.8  lb  at  60'  1.8  lb  at  800' 
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OOO-ft  source,  Weston's  values  have  been  plotted  in  Figure  21,  along  with  the 
three  Squaw  shots  and  the  Idealized  800-foot  source  level  curves  of  the  present 
,-tu.iy.  This  figure  indicates  that  whenever  the  two  quasi-theoretical  source  level 
:urvts  differ  significantly,  Weston's  analytical  curve  shows  less  satisfactory 
n&jceaent  vitr.  the  measured  data.  This,  together  with  the  generally  satisfactory 
comparisons  with  available  data,  is  some  indication  that  the  source  levels  of  this 
study  represent  an  Improvement  over  Weston's  analytical  formulation. 

8.  SIMHARY  AND  COKCUUSIONS 

8.1  The  computational  methods  developed  in  this  study  can  be  used  to  estimate 
acoustic  source  levels  for  small  charges  (weighing  up  to  a  fev  lb)  detonated  at 
depths  down  to  about  1000  ft.  Results  for  very  shallow  sources  (e.g.  60  ft)  will 
remain  less  reliable  than  those  for  deeper  sources  until  we  acquire  suitable  data 
to  compare  with  the  idealized  shallow-hurst  model.  The  quasl-theoretlcal  pressure- 
time  histories  and  their  corresponding  energy  levels  reported  here  are  believed  to 
be  the  best  estimates  available  at  this  time. 

8.2  Detailed  plots  (analysis  bandwidth  approximately  1  Hr)  of  idealized  source 
spectra  at  100  yards  range  are  shown  in  Figures  5-9  for  the  following  configurations: 
3-lb  charges  at  60  and  500-ft  detonation  depths,  and  1.8-lb  charges  at  60,  300,  and 
800- ft  detonation  depths.  Sstlmated  source  levels  in  l/3- octave  and  1-octave 
bandwldths  are  sunnarlzed  in  Tables  3  and  *♦,  respectively.  The  1.8-lb  results 
apply  to  underwater  sound  signals  (SUS  charges)  Mk  57,  Mk  61,  and  Mk  82  with 
appropriate  depth  settings;  the  3“lh  charge  results  correspond  to  test  conditions 
used  in  an  extensive  research  study. 

8.3  In  some  instances  considerable  difference  was  found  between  the 
theoretical  source  levels  presented  here  and  previous  measurements.  For  example, 
although  overall  agreement  with  the  data  of  Christian  (Refs.  2  and  25)  was  good, 
deviations  of  up  to  3  db  occurred  in  certain  frequency  bands  (paragraph  7.3,  7*5, 
7.6).  Significant  differences  from  Christian's  data  can  be  resolved,  since  all 
necessary  information  about  methods  used  in  data  recording  and  processing  is 
available  at  this  Laboratory.  Unfortunately,  the  same  la  not  true  for  the  data  of 
Turner  and  Scrlmger  (Ref.  12).  As  shown  in  Table  7,  the  theoretical  spsctnmi 
levels  differ  significantly  from  the  data  of  Ref.  12  for  low  frequencies- -almost 

5  db  in  two  cases .  These  deviations  may  be  real,  or  they  may  be  due  largely  to 
data  recording  and  processing  differences  idiich  can  easily  introduce  variations  of 
several  db  tlien  the  analysis  bandwidth  Is  narrow  relative  to  the  spectral  pattern. 

8.U  In  some  instances,  wide  discrepancies  (up  to  —  8  db)  exist  between  the 
source  level  estimates  of  the  present  study  and  calculations  using  the  analytical 
model  of  Weston1.  It  has  been  shown  that  an  inadequate  formula  for  the  shock  wave 
impulse  was  used  in  Weston's  low  frequency  model,  and  a  less  realistic  model  of  the 
explosion  wave  form  than  that  of  the  present  study  was  used  at  higher  frequencies. 
Limited  comparisons  with  available  data  indicate  the  present  computational  method 
yields  more  accurate  source  level  estimate;,  than  Weston's  analytical  model  *or  the 
800- ft  source. 

8.5  The  analysis  bandwidth  used  is  an  important  parameter  in  source  level 
definition.  Because  of  the  o?cillatory  character  of  explosion  spectra,  analyses 
in  1  Hz,  l/3-octave  and  octave  bands  may  yield  far  dlffer-r.t  source  level  values 
at  the  same  center  frequency. 
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8 .C  ceemiagly  snail  variations  in  burst  depth  can  cause  significant 
variations  in  narrow-band  source  level.  For  example,  experimental  variation  about 
tiie  300- ft  nominal  depth  setting  for  a  1.8- lb  SUS  may  produce  at  much  as  6.3  db 
source  M>el  difference  (in  tne  sO  nz  1/3-octave  band)  for  depth  variation*  of 
700  to  900  ft.  To  avoid  this  discrepancy,  the  depth  of  an  explosive  a>urce  must 
be  accurately  determined,  and  a  source  level  for  that  depth  determined.  No  simple 
correction  for  depth  variation  can  generally  be  applied  to  the  source  ltvel  to 
account  for  depth  variation. 

3.7  The  effect  of  the  surface  reflection,  sometimes  erroneously  included  as 
part  of  the  source  level,  has  been  demonstrated  in  one  instance  (paragraph  C.  10> 
to  produce  errors  in  the  source  level  at  different  frequencies  of  aa  much  ea  +U  and 
-5  db.  Since,  for  shallow  sources,  the  surface  reflection  cannot  be  separated 
from  the  direct  arrival,  direct  measurement  of  these  source  levels  is  not  possible. 
(The  source  level  is  not  necessarily  identical  to  the  measured  energy  spectrum. 

The  latter  nay  Include  surface  reflections  and  other  local  effects;  the  former 
3hculd  not. )  Furthermore,  no  simple  correction  of  a  constant  number  of  db  will 
correct  the  measured  spectrum  to  the  correct  source  level  values  at  all  frequencies. 

9.  RECOMMENDATIONS  FOR  FURTHER  WORK 


9*1  Further  comparisons  of  the  results  of  this  and  other  calculations  with 
shallow  data  must  be  made.  These  comparisons  would  evaluate  the  extrapolation  of 
shock  wave  parameters  at  the  shallow  depths  used  In  this  study,  as  veil  as  determine 
which  approach  is  best  suited  for  the  estimation  of  shallow  charge  source  levels. 
Suitable  data  would  consist  of  broadband  (~  0  -  20K  Hr.)  pressure-time  records, 
recorded  at  close  range  (preferrably  near  100  ysurdo).  These  would  allow  tha 
details  of  the  surface  reflection  to  be  analyzed  and  included  In  the  calculations 
for  comparisons  (similar  to  the  procedure  of  paragraph  7.6).  Narrow  band 
recordings  make  the  analysis  of  time  domain  characteristics  difficult.  Such  shallow 
data  may  be  available  at  other  facilities.  If  not,  suitable  field  work  should  be 
planned  to  gather  the  needed  data. 

9.2  In  order  to  evaluate  the  effect  of  different  analysis  techniques,  aeverml 
samples  of  data  should  be  processed  through  different  analysis  systems,  digital, 
analog  and  hybrid  digital-analog,  and  the  results  compared.  Since  our  analysis 
capability  Is  digital,  we  would  be  eager  to  exchange  data  with  facilities  using 
other  techniques. 

9.3  la  the  absence  of  the  above  additional  comparisons,  we  tentatively 
recommend  that  shallow  explosive  source  levels,  idiere  the  surface  reflection  cannot 
be  separated  out,  be  determined  by  the  method  of  the  present  study.  As  direct 
measurement  of  these  source  levels  is  impossible,  we  feel  that  the  promising 
comparisons  of  section  7  indicate  that  this  method  yields  the  most  accurate 
available  estimates  for  such  shallow  sources.  Close-in  pressure  measurements  are 
still  needed  In  particular  applications  so  that  the  actual  source  depth  may  be 
determined,  and  the  error  introduced  by  source  depth  variation  eliminated. 
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FIG.  6.  IDEALIZED  SOURCE  LEVEL-3  LB  OF  TNT  AT  500  FOOT  DEPTH 
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IDEALIZED  SOURCE  LEVEL-1.8  LB  OF  TNT  AT  300  FOOT  DEPTH 


OF  TNT  AT  800  FOOT  DEPTH 


OCTAVE  AND  l/3-OCTAVE  BAND  SOURCE  LEVELS- 3  POUNDS 
OF  TNT  AT  60  FOOT  BURST  DEPTH 
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TNT  AT  300  FOOT  BURST  DEPTH 


OF  TNT  AT  800  FEET  BURST  DEPTH 


FOOT  DEPTH  WITH  SURFACE  REFLECTION  AT  24  MSEC 


16.  DETAILED  ENERGY  SPECTRA  SHOWING  EFFECT  OF  SURFACE  REFLECTION 
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FIG.  20.  OCTAVE  BAND  COMPARISON  WITH  CHRISTIAN’S  SOURCE  LEVELS  -  800  FOOT  DEPTH 


ENERGY  FLUX  DENSITY 
(dB  re  1  erg/cm^/Hz) 
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- 1.8  LB  AT  80C  FT-IDEALIZED  SOURCE  LEVEL 

X - X  WESTON'S  ANALYTICAL  FORMULATION 


A  SQUAW  SHOT  222-2-1 0 
O  SQAW  SHOT  323-2-10 
•  SQAW  SHOT  76-3-10 


FREQUENCY,  Hz 

FIG.  21.  ONE-THIRD  OCTAVE  BAND  SOURCE  LEVEL  COMPARISONS  -  800  FOOT  BURST  DEPTH 
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IMPULSE  FROM  REF.  16  =  |$ 
IMPULSE  FROM  REF.  18  =  I  w 


FREQUENCY,  Hz. 

COMPARISON  WITH  WESTON'S  ANALYTICAL  FORMULATION 
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